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We study the effect of buckling and breathing phonons, relevant for cuprate superconductors, 
on the d-wave superconductivity in the two-dimensional Hubbard model by employing dynamical 
cluster Monte Carlo calculations. The interplay of electronic correlations and the electron-phonon 
interaction produces two competing effects, an enhancement of the effective d-wave pairing inter- 
action, which favors d-wave superconductivity, and a strong renormalization of the single-particle 
propagator, which suppress superconductivity. Due to the later effect we find that buckling and 
breathing phonons suppress the superconductivity in the region of parameter space relevant for 
cuprate superconductors. 



Introduction. Many experiments, includingRaman 
neutron scattering and photoemission [3,^ find clear 
evidence of electron-phonon (EP) interaction in high Tc 
materials. While it is widely accepted that strong elec- 
tronic correlations play a fundamental role in the mech- 
anism of high temperature superconductivity, the role of 
EP interaction is under hot debate. 

The phonons believed to be the most relevant to the 
cuprates are the buckling (BC) and the breathing (BR) 
modes [H, [3, 0, IE S IJi defined, respectively, by the 
out-of-plane and in-plane displacement of the oxygen 
ions. Previous investigations suggest that the symmetry 
of these phonon modes is significant. There are argu- 
ments suggesting that coupling to the out-of-plane oxy- 
gen BC or the local Holstein (H) phonons, enhance the 
d-wave pairing interaction [1, [3, Elj IiSIj whereas, 
the BR mode has been found to suppress the interac- 
tion 

B E, El- It also was shown that the electronic 
correlations in the presence of EP coupling strongly en- 
hance polaron formation l^, 3, l3| and renormalize the 
quasiparticle (QP) weight, an effect which suppresses su- 
perconductivity. 

In this letter we employ the dynamical cluster approx- 
imation (DC A) 1^, 17 1 to study two models relevant for 
cuprates, the Hubbard model with BC and BR phonons. 
We also extend the work on superconductivity in the 
Hubbard- Holstein model [l^. DCA has proved to be 
one of the few techniques able to capture the supercon- 
ducting properties of strongly correlated systems 17 1. 



We investigate the role of phonons on superconductiv- 
ity in the region of small and intermediate doping where 
the antiferromagnetic (AF) correlations are strong. In 
agreement with previous investigations, we find enhance- 
ment of the d-wave pairing interaction by H [1, 15| and 
BC [i. Hoi , [ill . IT^ phonons. Moreover we find an en- 
hancement of the d-wave pairing interaction even for BR 
phonons contrary to some predictions [13, U | . How- 



ever, despite the enhancement of the d-wave pairing, we 
find a suppression of the superconducting Tc for all three 
modes, due to the strong renormalization of the elec- 
tronic single-particle propagator. 

Formalism. The Hamiltonian for each model, 
Hubbard-Holstein (HH), Hubbard-Buckhng (HBC) and 
Hubbard-Breathing (HBR), can be written as 



H = Hr 
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is the Hubbard part with nearest-neighbor hopping t and 
on-site repulsion U . For H phonons, described as a set 
of independent oscillators at every site i which couple 
locally to the electronic density. 
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where {ut, Pi} are canonical conjugate coordinates for 
each oscillator with characteristic frequency luq . The BC 
(BR) phonons are independent out-of-plane (in-plane) os- 
cillators on each bond i + 7/2, with j ^ x,y, such that 
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where the last term of Eq. |4] with plus (minus) describes 
the coupling to the BC (BR) phonons. The dimen- 
sionless EP coupling, defined as the ratio of the single- 
electron lattice deformation energy to half of the elec- 
tronic bandwidth W/2 = At, is = 2g'^ / {2MujIW) and 
\BC,BR ^ 4 X 2g^/{2MujlW) [3], with an extra factor 
of four in X^c.br local coordination. Note that in 
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general the coupling to the BC mode implies both mod- 
ulation of orbital energy and hopping integrals [3| ■ How- 
ever, due to the technical difficulties associated with off 
diagonal EP coupling, we consider here only the former 
effect as was done previously in Ref . 0, [13] ■ 

To study the Hamiltonian ^ we employ the DCA, 
a cluster mean-field theory which for a two dimensional 
system maps the original lattice model onto a periodic 
cluster of size Nc = Lj. embedded in a self-consistent 
host. Correlations up to a range Lc are treated explic- 
itly, while those at longer length scales are described at 
the mean-field level. With increasing cluster size, the 
DCA systematically interpolates between the single-site 
dynamical mean field result and the exact result, while 
remaining in the thermodynamic limit. Cluster mean 
field calculations on the simple Hubbard model success- 
fully reproduce many of the features of the cuprates, in- 
cluding a Mott gap and strong AF correlations, d-wave 
superconductivity and pseudogap behavior [l3| . 

We solve the cluster problem using a quantum Monte 
Carlo (QMC) algorithm modified to perform the sum 
over both the discrete field used to decouple the Hubbard 
repulsion, as well as the phonon field u. The space of con- 
figurations of the latter field is significantly larger than 
the former. In part, this is offset by the strong correla- 
tions of the phonon field in Matsubara time. Therefore, 
in the QMC Markov process, correlated changes, in which 
adjacent phonon fields in time are moved together, are 
mixed in with local moves to reduce the autocorrelation 
time. Nevertheless, the present code including the effect 
of phonons requires significantly more CPU time than re- 
quired for the Hubbard model. Thus, most of the present 
calculations are restricted to clusters of size iVc = 4. Nev- 
ertheless we check the robustness of our conclusions by 
employing calculations on larger clusters of size Nc = 16. 

Results. In a previous paper, addressing the HH 
model [3] , we show that the synergistic interplay of AF 
correlations and the EP interaction strongly enhances 
both polaron formation and antiferromagnetism. Since 
the effective EP coupling is inversely proportional to the 
kinetic energy of the holes, the AF correlations, which 
reduce the mobility of the holes, enhance the effective 
EP coupling. On the other hand, since at finite doping 
the antiferromagnetism is suppressed by the hole mo- 
tion, the decrease in hole mobility due to EP coupling 
increases the antiferromagnetism. Considering the BC 
and the BR phonons we find a similar behavior indica- 
tive of enhanced polaron formation. In Fig.[T]-a we show 
that the AF susceptibility at 5% doping is strongly en- 
hanced by coupling with any of the three phonon modes. 
The EP coupling also reduces the Matsubara QP weight, 
Zo(T) = 1/(1 - lmT,{K,iTTT)/7:T) (Fig. [1] -b for mo- 
mentum K — (0,7r)), and increases the kinetic energy 
(Fig. [T]-c) [2^. For all three models the local moment, 
/i^ = ((n| — at low temperature has a large value, 

close to that corresponding to the A = case (Fig. [T] 
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FIG. 1: (color online) a) AF susceptibility xaf, b) Matsubara 
QP fraction Zq = 1/(1 - ImY;{K,TTT)/-KT) at K = (0,7r), c) 
kinetic energy Ekin = X/fco- ^{^)^k,cr and d) local moment 
/^^ = ((n| - Uif) versus T for HH, HBC and HBR models at 
5% doping. All three modes show enhanced polaron formation 
at low T. 
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FIG. 2: (color online) Tc versus A for HH, HBC and HBR 
models at (a) 5% and (b) 15% doping for ujo — 0.3t, U — 8t 
and Nc = 4. 



-d), due to reduction of the effective hole hopping. How- 
ever, aside from the similarities between the three mod- 
els, there are also significant differences. We find that for 
the HBC model the same value of EP coupling produce 
much stronger effects then for the HH and HBR models. 
For example at 5% doping a value of A^'-^ = 0.3 produces 
effects similar to X^'^^ = 0.5. The temperature depen- 
dence of the local moment (Fig. [T]-d) as well as of other 
quantities such as DOS (not shown) is much stronger in 
the HBR model. 

In Fig. [2] we show the dependence of superconducting 
Tc versus EP coupling for the three models for ujq — 
0.3t at 5% and 15% doping. We find Tc decreases with 
increasing A, weakly for small A and then quite abru ptly 
for A > Ac. We find that \^,bc,br ^ 0.5,0.25,0.4 
Notice that the decrease of Tr with A for HBC model is 
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much sharper once A > Af*-^. 

Previous investigations predict an enhancement of the 
phonon contribution to the d-wave pairing interaction 
when the coupUng with BC or H phonons is consid- 
ered [1, Q- However, an increase of the interaction which 
favors d-wave pairing does not necessarily imply an in- 
crease of Tc, since the reduction of the QP weight and 
DOS at the Fermi level, N{0), has the opposite effect. 
We find for all three models an increase of the d-wave 
pairing interaction. However, in spite of this, the Tc is re- 
duced as shown in Fig. ^ To understand this dichotomy, 
we must disentangle the effects of the pairing interaction 
from renormalization of the single-particle propagator. 

A divergent pairing susceptibility indicates the super- 
conducting transition and hence Tc- At this instabil- 
ity the leading eigenvalue of the pairing matrix becomes 
1. The pairing matrix M = Txo enters in the Bethe- 
Salpeter equation of the two-particle Green's function. 
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X = Xo + XoTx = Xo + Xo(l - M)-H4 . 
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where xo = G * G (the bubble diagram) describes prop- 
agation of the two particles without mutual interac- 
tion, a product of the fully renormalized single-particle 
propagators G{k,iuu) — [iuo — e(fc) -I- ^ — Yj{iLu,k))~^ . 
^Q=o,iv=o{.K,iuj]K' ,iuj') is the irreducible interaction 
vertex in the particle-particle channel and can be re- 
garded as the effective renormalized interaction. We 
find that for all three models the leading eigenvector 
i^) of the pairing matrix M has d-wave symmetry. 
The phonons affect the pairing matrix by modifying 
both the effective interaction F and the single particle 
propagator G, or implicitly xo- To separate these ef- 
fects we compare, in FigI3J the eigenvalues of two dif- 
ferent pairing matrices: M and Mq. M was defined 
above and we define Mq = Fxoo with xoo = G° * G°. 
G^{k,iuj) — {iuj — e(fc) ^ \iu — T,u{iuj, k))~-^ is a single- 
particle propagator which does not account for the renor- 
malization resulting from the EP interaction, i.e fijj and 
S(7 are the obtained by setting A = 0. Note that the ef- 
fective interaction F in Mq is, however, fully renormalized 
by phonons. 

The leading eigenvalue of M increases with decreasing 
T and reaches a value of 1 at the transition temperature 
Tc- This transition occurs at lower temperature when A 
is finite, thus indicating a reduction in Tc in accordance 
with the phase diagram shown in Fig. [2l However, the 
d-wave eigenvalue of the matrix Mq increases much faster 
and reaches 1 at a temperature much larger than without 
phonons. Since the matrix Mq contains the fully renor- 
malized interaction and the single-particle propagators 
which are not renormalized by phonons, this shows that 
the EP coupling strongly enhances the d-wave pairing 
interaction. However, the competing effect, renormaliza- 
tion of XO: is also very strong such that the net effect is 
a reduction of Tc. 
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FIG. 3: (color online) d-wave eigenvalue of pairing matrices 
M = r xo and Mo = F xoo versus T on iVc = 4 cluster for HH 
(a), HBC (b) and HBR (c) models at 5% doping and HBC (d) 
at 15% doping. Despite the overall reduction of Tc shown by 
the eigenvalues of AI the eigenvalues of A/q increase strongly 
with decreasing T, showing that EP coupling enhances the 
effective pairing interaction. 



Alternatively, to investigate the pairing in- 
teraction one can look at quantities such as 

and PdQ = J2K.iuj^diKi'^'^)Xo{K,iuj) which are the 
respective projections of the interaction vertex and xo 
on the subspace spanned by the d-wave eigenvector. 
These were previously defined in Ref . [S^l . Fig. 3] -a and 
-b show that EP coupling enhances the effective pairing 
interaction Vd and strongly reduces P^o, reinforcing the 
conclusions previously drawn from Fig. [3l 

Note that we see an enhancement of d-wave pairing in- 
teraction for all three phonon modes, including BR. DCA 
calculations in the Hubbard model without phonons show 
that the main contribution to the d-wave pairing is 
contained in the particle-hole spin S — 1 channel at 
Q — (tt, tt) (21| , i.e. the pairing is a result of exchanging 
AF spin fluctuations. We speculate that the increase of 
the pairing interaction when EP is present results from 
the enhancement of the AF susceptibility. However a 
decomposition of the pairing vertex F in the fully ir- 
reducible and partially reducible particle-hole spin and 



density components, similar to the one done in Ref. [21 1 
for the Hubbard model, is necessary to better understand 
the effect of phonons on the pairing interaction. 

It is important to ask whether these results will change 
for other regions of parameter space relevant for cuprates 
or when larger clusters are considered. We find that EP 
coupling always reduces Tc for small phonon frequency 
u!q < t. We exemplify this by showing in Fig|3]-c, -d and 
-e the inverse of the d-wave pairing susceptibility Xd for 
some particular cases. Larger cluster, Nc = 16 sites, and 
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rithm we investigate the Hubbard model with H, BC 
and BR phonons. We find that the interplay of EP in- 
teraction and electronic correlations in HBC and HBR 
models leads to synergistic enhancement of both polaron 
formation and AF correlations, similar to previous find- 
ings for the HH model ■ Regarding superconductivity, 
all three modes produce two competing effects: a strong 
renormalization of the single-particle electronic propaga- 
tor which suppresses superconductivity and an enhance- 
ment of the effective pairing interaction which favors d- 
wave superconductivity. In the region of parameter space 
relevant for cuprates, we find that the combination of 
these two effects leads to a reduction in superconducting 



FIG. 4: (color online) a) d-wave pairing interaction Vd and 
b) d-wave projected bare bubble Pdo versus T for U = 8t, 5% 
doping, Luo ~ 0.3t and A'^c = 4 cluster. Vd (Pdo) is enhanced 
(reduced) by EP coupling, c) Inverse of the d-wave pairing 
susceptibility 1/Xd for HH, HBC and HBR models for U — 8t, 
Nc = 16 cluster and luq = 0.15t at 15% doping, d) 1/xd 
for HBC model for U = 8t, loq = 0.3t, Nc = 4 cluster at 
25% doping, e) 1/xd for HBC model for U = 4t, loq = Q2t, 
Nc = 4 cluster at 5% doping. We find that a finite A always 
suppresses Tc- 



smaller phonon frequency, ujo = 0.15t, results are illus- 
trated in FigU-c for all three models. Even though we 
were unable to reach temperatures equal to Tc, due to 
the sign problem present in large clusters calculations, it 
is obvious that a finite A suppresses the pairing suscep- 
tibility. Regarding the phonon frequency dependence of 
Tc, we find that, in general, Tc decreases with decreasing 
luq (not shown), i.e. a positive isotope effect. However 
the magnitude of the isotope effect depends strongly on 
the mode, A and doping. A detailed investigation will 
be presented in a subsequent publication. In Fig. 2] -d 
and -e we show that EP coupling reduces Tc for the HBC 
model in the overdoped region, 25% doping, and in the 
weak coupling regime, U — At. The same conclusion can 
be drawn for the HH and HBR models (not shown). 

Our results indicate that the coupling of the electronic 
density to local phonons is, in general, not favorable for 
superconductivity in the cuprates. This result is contrary 
to some previous speculation [1, [i^ about the role played 
by the EP interaction in high Tc, and therefore signifi- 
cant to the experimental community. Synthesis aimed at 
increasing Tc by designing new materials with strong EP 
coupling may lead to the antithetical result. A proper 
treatment of quasiparticle dressing is key to understand- 
ing the role played by the EP interaction. Methods that 
fail to properly address this point, e.g. by truncating the 
phonon Hilbert space, may lead to incorrect or erroneous 
conclusions. 

Conclusions. By employing DCA with a QMC algo- 
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